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ABSTRACT 

Our Sun, primarily composed of ionized hydrogen and helium, has a surface temper¬ 
ature of 5777 K and a radius Rq ss 696,000 km. In the outer Rq/3, energy transport is 
accomplished primarily by convection. Using typical convective velocities u ~ 100 ms“^ 
and kinematic viscosities of order 10“^ m^s“^, we obtain a Reynolds number Re ~ 10^^. 
Convection is thus turbulent, causing a vast range of scales to be excited. The Prandtl 
number, Pr, of the convecting fluid is very low, of order 10“'^ -10“^, so that the Rayleigh 
number (~ Re^Pr) is on the order of 10^^ — 10^^. Solar convection thus lies in ex¬ 
traordinary regime of dynamical parameters, highly untypical of fluid flows on Earth. 
Convective processes in the Sun drive global fluid circulations and magnetic fields, which 
in turn affect its visible outer layers (“solar activity”) and, more broadly, the heliosphere 
(“space weather”). The precise determination of the depth of solar convection zone, de¬ 
partures from adiabaticity of the temperature gradient, and the internal rotation rate as 
a function of latitude and depth are among the seminal contributions of helioseismology 
towards understanding convection in the Sun. Contemporary helioseismology, which is 
focused on inferring the properties of three-dimensional convective features, suggests 
that transport velocities are substantially smaller than theoretical predictions. Fur¬ 
thermore, helioseismology provides important constraints on the anisotropic Reynolds 
stresses that control the global dynamics of the solar convection zone. This review 
discusses the state of our understanding of convection in the Sun, with a focus on helio- 
seismic diagnostics. We present our considerations with the interests of fluid dynamicists 
in mind. 

Subjeet headings: Helioseismology, Convection, Non-Boussinesq, Dynamo 


1. Introduction 


Solar energy is powered by nuclear fusion in the Sun’s core (inner 20% by radius) through p-p 
chain fusion, whereby ionized hydrogen is converted to helium. Other by-products of fusion include 
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high-energy photons and neutrinos. Density and temperature are very high in the core, on the order 
of 150 g/cvo? and 15.7 million K, respectively. High-energy photons released in fusion reactions 
are absorbed in a few millimeters of solar plasma and then re-emitted in random directions and at 
lower energies {free-free scattering). Photons gradually diffuse outwards from the core to the Sun’s 
surface; random walk estimates suggest a transit time in excess of 100,000 years (while neutrinos 
take only a few seconds to reach the surface). Since the temperature gradient is sufficiently small 
below about O.7i?0, radiative transfer by photons is the primary means of energy transport from 
the core. At the base of the convection region, i.e. r ~ 0.7Rq, the temperature is approximately 2 
million K. Radiation becomes less efficient further outwards because of the increase in the opacity 
of the plasma, and much of the energy transport occurs by thermal convection. The radiative 
interior and the convective zone are separated by a thin layer (~ O.O5i?0, e.g. Elliott &: Gough 


1999; Christensen-Dalsgaard et al. 2011) where the stratification changes rapidly from convective 


stability to marginal instability. This region also shows a relatively sharp change between the solid- 
body rotation of the radiative interior and the differential rotation of the exterior convection zone; 
it is termed the tachocline ( Spiegel Sz Zahn||1992 ). 

Convective motions are thought to drive the dynamics and activity of the Sun. On small scales, 
convection directly generates and sustains magnetic fields (the Sun’s small-scale dynamo field; 
e.g. Vogler et al. 2005). On larger scales, the turbulent convection in the rotating Sun produces 


Reynolds stresses which drive large-scale (global) motions. These large-scale motions amplify the 
global magnetic field whose manifestations at the solar surface are sunspots of magnetic complexes, 
which in turn are stressed by convective motions, producing flares and space weather. Magnetic 
fields act on long timescales to mediate angular momentum loss from the Sun as a whole (magnetic 
braking). Convection also mixes plasma inside the Sun; this type of mixing can significantly affect 
the lifetimes of main-sequence stars (that is, stars that are fusing hydrogen atoms to form helium 
atoms in their cores). For these reasons it is important to understand convection inside stars, 
particularly the Sun. 


Convective motions in the Sun are more complex than in laboratory experiments of Rayleigh- 
Benard convection, in particular because of stratification. For instance, the density contrast be¬ 
tween the base and the top of the solar convection zone is about 10®. Furthermore, the energy 
transport by convection is very efficient in the Sun, in that only very small departures from the 
adiabatic temperature gradient are required. These motions are also difficult to model in faithful 
detail: the relevant parameters in the Sun are extreme, with the highly turbulent flows involving 
spatial and temporal scales spanning many orders in magnitude. For reviews on solar convec¬ 
tion, see e.g. Miesch (2005), Miesch Sz Toomre (2009), Rieutord et al. (2010), and Hanasoge & 


Sreenivasan (2014). 


In this context observational tests of our understanding are essential, and we shall focus in this 
review on the contributions made by helioseismology. The accurate inference of the depth of the 
convection zone, the measurement of the departure from adiabaticity as a function of depth, and the 
determination of solar interior rotation and meridional circulation (at least in the surface layers) are 
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all important contributions of tielioseismology (Ctiristensen-Dalsgaard 2002; Gizon & Birch 2005). 
The current focus is on characterizing the three-dimensional (3-D) properties of convective flows 


and magnetism in the solar interior using techniques of local helioseismology (Gizon, Birch &: Spruit 


2010 ). 


1.1. Comparison with laboratory convection 


Standard experimental work on convection is carried out in the Rayleigh-Benard apparatus. 
Nominally, a shallow fluid layer contained in a cylinder with non-conducting side walls and con¬ 
ducting bottom and top walls is heated from below and cooled from above, so that the entire heat 
from the bottom wall goes through the fluid to the top wall. An important question concerns the 
amount of heat that is transported across the fluid for a given temperature difference maintained 
between bottom and top walls. Heat is transported purely by molecular conduction if the temper¬ 
ature difference is small (in a sense to be clarihed), but the fluid begins to move and enhance the 
heat transport thereafter. A non-dimensional measure of the heat transport is the Nusselt number, 
which is the ratio of the actual heat transport to that possible by conduction alone. The buoyancy 
effect produced by the temperature difference is balanced by the smearing effects of diffusion and 
viscosity, the ratio of the former to the latter being called the Rayleigh number; the larger the 
Rayleigh number the larger is the effective temperature difference. One goal of convection studies 
is to discover how the Nusselt number scales with the Rayleigh number. Much is known about this 
flow ( Gastaing et al.||1989 Niemela &: Sreenivasan||2006b Ahlers||2009 Ghilla fc: Schumacher||2012 ) 


though the nature of the heat transport law for asymptotically large Rayleigh numbers is still a 


topic of lively discussion ( 

Siggia||1994 

Niemela Sz Sreenivasan||2003 

2010 

Ahlers et al.||2012 

Urban 

et al. 

2014; 

Stevens, Glercx & Lohse 

2013 

Shishkina et al. 

2015 

)• 


It is also possible to establish a flow in which one fixes the rate of heat transport at the bottom 
and top walls, instead of their temperatures. It is this latter condition that might apply to the 
Sun and other stars. The convection boundaries are evidently not solid but steep changes in the 
adiabatic temperature gradient that occur at the base of convection zone signihcantly restrict free 
movement of plasma across this layer. There can be small-scale penetrative convection; indeed, 
it is thought to play a critical role in stellar rotation spin down and, therefore, measuring the 


depth of penetration is an active area of research in helio- and astero-seismology {e.g. Ghristensen- 


Dalsgaard et al. 2011). Nevertheless, all large scale motions are restricted at the base resulting 


in a plausible analogy with a solid surface. Likewise, at the top of the convection region, steep 
gradients in temperature and density caused within the near-surface convection region endow it 
with boundary-like properties. 


Thermal energy produced in the core and transmitted at this radiative-convective interface 
provides the constant heat flux boundary condition at the bottom of convection zone. However, 
the variation of thermodynamic properties and the gas composition between the top and bottom 
of the convection zone renders stellar convection highly non-Boussinesq. We recall that only small 


























































- 4 - 


deviations from the large adiabatic gradients of temperature and density matter for convection. 
Yet the Rayleigh numbers are vast and motions inside the convection zone are turbulent. 


It is useful to note that the Prandtl number of the fluid is very small everywhere in the con¬ 
vective flow regime of the Sun. Even fluids like mercury and liquid sodium, sometimes used in the 
laboratory, possess Prandtl numbers which are several orders of magnitude larger than solar plasma. 
This could have important consequences, but the role of the Prandtl number at very high Rayleigh 
numbers is not well understood. For small Prandtl-number fluids, no significant heat transfer en- 


hancement is seen in laboratory measurements ( 

Ecke & Niemela 2014 

Ahlers 

2009 

) and direct 

nnmerical simulations ( 

Horn &: Shishkina|2015 

Stevens, Glercx &: Lohse||2010 

Oresta, Stringano & 


Verzicco 2007) of rotating turbulent Rayleigh-Benard convection. Ekman vortices develop irregu¬ 


larly and are shorter than in the case of large Prandtl-number fluids. The thermal diffusivity, larger 


than the kinematic viscosity, reduces the effectiveness of Ekman pumping (Stevens, Clercx & Lohse 


2010; Horn & Shishkina 2015). For all Prandtl numbers, Horn &: Shishkina (2015) identify the 


different regimes of rotating Rayleigh-Benard convection by comparing the poloidal and toroidal 
energies of the flow field. This identification was found to be robust to changes in the aspect ratio 
and geometry of the container. 


While the entire region of the Sun beyond r ~ 0.7i?Q is convective, it may be conceptually 
useful here to consider the outer 0.3Rq as two separate regions. “Near-surface convection” is a 
thermal boundary layer with a thickness on the order of 300 km, where radiative opacities change 
rapidly and the plasma transitions from optically thick to thin, i.e., it becomes more transparent. 
The rest of the convection region, forming the bulk, is termed the “deep convection” region. These 
two regions of the convection zone are likely connected by the cool descending plumes associated 
with overturning surface convection, seeded at the photosphere. The fate of these plumes, i.e. their 
overall stability during descent sets the degree of mixing and thus the dynamics of global convection 
{e.g. 


Spruit||1997 Nordlund, Stein & Asplund 2009 


1.2. Near-surface convection 


Near-surface convection, marked by the presence of convective cells of characteristic size 1000- 
2000 km called granules, can be observed by high-spatial-resolution optical imaging devices (see 


Nordlund, Stein & Asplund 2009, and references therein). The centres of granules are hotter than 
the edges, implying that granulation is a form of overturning convection, where hot upward moving 
plasma radiates away the thermal energy, and channels between granules (intergranular lanes) 
provide conduits for the cool fluid to reenter the solar interior. Granules are the primary vehicles of 
thermal transport in the near-surface layers of the Sun, carrying nearly the entire solar luminosity 
outwards ( Nordlund, Stein Sz Asplund]2009 ). Granules have been the subject of careful observation 
{e.g. Schrijver, Hagenaar &: Title||1997) and properties such as their size and intensity-contrast are 


well characterized. Gomprehensive numerical simulations of surface convection (Stein &: Nordlund 


2000; Vogler et al. 2005) are successful in reproducing these properties. 
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Fig. 1.— Line-of-sight projected velocity field measured at the photosphere of the Sun by the 
Michelson Doppler Imager. Average rotation has been subtracted. The structure visible in this 
image is supergranulation, a flow system of size 35000-km in which horizontal velocities are about 
300 m/s and radial velocities are ~ 30 m/s at the surface. Courtesy M. Rieutord, see Rieutord 


et al. (2010). 









- 6 - 


One challenge of the simulations concerns supergranulation, a distinct pattern of horizontal 
flows (root-mean-square velocity on the order of 300 m/s) with a characteristic size of 35 000 km 
and a lifetime of 1 ~2 days, which is easily seen away from the disk center in direct Doppler images 
(see Figure [T| |Rieutord &: Rincon 2010a; Hathaway et al. 2000). Although supergranulation was 
discovered by [Hart (1954) more than sixty years ago, its origin and subsurface properties still 
pose major challenges. It has been suggested that it may correspond to a preferred scale of deep 


convection 

Leighton, Noyes &: Simon 

1962 

) or to a nonlinear dynamical interaction of granules 

(Hast 

2003 

Rieutord &: Rincon 

2010b 

. The supergranulation pattern propagates faster than the 


solar plasma and has been cited as an example of traveling-wave convection (Gizon, Duvall & Schou 


2003). 


1.3. Deep convection 


Convection in the interior is a difficult problem to address because it is both optically inac¬ 
cessible and the governing physics lies in a parameter regime that is beyond laboratory studies or 
numerical simulations. The basic issue of whether deep convection is controlled by heating from 


below or by the radiative cooling from above or by a combination of the two is not clear (Spruit 


1997). Owing to steep stratification (nonetheless far smaller than that in the near-surface convec¬ 
tion region), fluid properties such as the kinematic viscosity and thermal and magnetic diffusivities 
at the base and top of the convection zone, and thus the time and length scales in those regions, 
are dramatically different. Therefore if heating and cooling were to be decoupled, i.e. if mixing 
of descending plumes into the surrounding medium were weak, it is conceivable that the problem 
would separate into more tractable sub-problems. Whereas, if the coupling were strong, the physics 
of the convection zone would be altogether stiff since both small and large scales would play crucial 
roles. The physics of the coupling of convection at the surface thermal boundary layer with deeper 
layers is thus less easy to understand than in standard Boussinesq convection. 

Numerical simulations typically show large-scale overturning convective cells, coherent on scales 
much greater than supergranulation. In other words, the surface convective power spectrum ob¬ 
tained from numerical calculations suggest that the largest scales are most energetic, in contrast to 


that on the surface. Simulations show elongated “banana cells” at low latitudes {e.g. Miesch et al 


2008), which have not (yet) been detected on the Sun. Using measurements of the photospheric flow 


field, Hathaway, Upton &: Colegrove (2013) measured large-scale fluid motions (giant cells), with 
velocities of the order of 8- 20 m/s, smaller by at least an order in magnitude than those associated 
with intermediate and other smaller scale convective structures (also see Figure]^. Giant cells are 
seen at all latitudes, with a seeming preference for high ones. 

One might speculate that the limited parameter regime of numerical calculations can only lead 
to standard turbulent mixing and conventional overturning convection. The fundamental question 
is whether, in the dramatically different regime prevalent in stellar interiors (that of the Sun in 


particular), the highly asymmetric up and downflows are well mixed at all (also see section 3.3) 
















































Differential rotation in the convective envelope provides a stringent test of simulations, and prior 


efforts have failed to reproduce it accurately from first principles {e.g. Balbus & Schaan 2012 


Miesch et al.|2012 ). Another constraint on the quality of the computed solutions is the prediction of 
meridional circulation, although uncertainties in our knowledge of meridional flow limit its potency 


as a constraint {e.g. Gizon, Birch &: Spruit 2010). In spite of inadequate resolutions, numerical 
simulations provide the best and detailed estimates of the properties of convection and dynamo 
action. 

The base of the convection zone contains a layer where the temperature gradient sharply transi¬ 
tions to strong stability. The motion of fluid is thus restricted and small-scale plumes penetrate into 
the overstable radiative interior. Convective penetration affects the mixing of chemical elements 
in the interior. It also excites internal gravity waves which in turn contribute to the transport 


of angular momentum and to the secular spin down of the Sun (Mathis 2013). The physics of 


convective overshoot is complicated by the presence of magnetism and the strong coupling between 
enormously different timescales in the convection zone and radiative interior. Appreciating the 


dynamics of this layer is actively pursued through the use of phenomenological models {e.g. Zahn 
1991), numerical simulations {e.g. Rogers & Glatzmaier 2005 Brun, Miesch & Toomre 2011), and 


seismic analyses {e.g. Christensen-Dalsgaard et al. 2011). 


2. Models of convection 
2.1. Mixing-length theory 


Since convection plays a crucial role in stellar evolution (by the mixing of elements), simplified 
prescriptions for its properties are of considerable importance. A widely used model is the mixing- 
length theory (MLT) of convection, originally due to Prandtl (1925) for hydrodynamic turbulence. 


extended to stellar physics by Bohm-Vitense ( 

1958 

)• 

the solar convection zone ( 

Spruit 

1974; 

Gougl 



codes. MLT is based on the simple phenomenology of parcels of fluid transporting thermal energy 
over a specihed length scale and subsequently immediately undergoing complete mixing. On the 
other hand, turbulence and convection are multi-scale phenomena. This crucial aspect of MLT, 
namely a single-scale description, does not handle the physics adequately {e.g. Weiss et al. 2004). 
Further, MLT does not model strong compressibility and interaction with radiation associated 
with near-surface layers, making it inaccurate there. Nevertheless, it provides a basic recipe for 
incorporating convection into stellar models and is therefore of great utility. More sophisticated 


modeling is required to address multi-scale dynamics and radiative transfer (for instance, see Ganuto 


&: Dubovikov 1998 for more detailed models of convection). 
























































Fig. 2.— Power spectrum of surface convection as a function of spherical-harmonic degree, as 
measured in line-of-sight Doppler velocity by the Helioseismic and Magnetic Imager. Granulation 
energy is broadly spread up to £ ~ 4000. A distinct peak corresponding to supergranulation is seen 
around i ~ 120. At low degrees, convective velocities fall sharply with decreasing i. Courtesy of 
David Hathaway. 
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Fig. 3.— Streamlines from a fully-compressible numerical simulation of near-surface convection 
(Stagger code; Nordlund, Stein &: Asplund 2009). The computational domain (a box 48 Mm on 
the side and 20 Mm deep) is seen from the bottom. On the left are the streamlines for up flows 
measured in the photosphere, on the right for downflows. The color along the streamlines indicates 
depth (white for deep, dark red for the surface). The downflows that reach the bottom of the box 
delineate the borders of supergranules. The upflows that originate from the deep layers are found 
at the center of supergranules. Courtesy of Robert Stein. 
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Fig. 4.— Temperature, radial vorticity, and horizontal divergence near the top (r = 0.98Rq) of a 
simulation of deep solar convection in a spherical shell (ASH simulation). The equator is indicated 
with a solid line. Magnified areas shown in the lower panels are square regions of size 45°. From 


Miesch et ah (2008). 
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2.2. Numerical simulations 


The surface layers of the Sun exhibit dynamic interplay between small-scale magnetism and 
granulation. Observations of granulation show a flow field that has distinct, seemingly well formed 
convective cells of typical size 1000-2000 km, with a lifetime of about 10 minutes. Fibril-like 


magnetic helds are swept to the intergranular lanes by horizontal flow helds {e.g. Pietarila Graham, 


Cameron & Schussler|2010 

). Granulation is also the primary source of acoustic oscillations observed 

on Sun’s surface {e.g. 

Goldreich & Keeley 

1977; 

Nordlund, Stein & Asplund 

2009 

). Simulating 


near-surface convection requires modelling radiative transfer and compressible hydrodynamics in a 
strongly stratified medium. Given these complexities, the accurate reproduction of the observables 
(such as granule size distribution, flow magnitude, and line formation properties) by simulations 
is remarkable (Stein &: Nordlund 2000| Vogler et ah 2005; Rempel, Schiissler & Knolker 2009[ 


Nordlund, Stein & Asplund |2009 e.g.)\ this is attributable perhaps to the steep density stratification 
(scale height of ~ 200 km), which causes the flow held to expand rapidly with height, smoothing 
small-scale huctuations and limiting the degree of interaction between the large and small scales. 

The computation of deep convection is a much more challenging problem, requiring one to 
incorporate a vast range of spatial and temporal scales, inaccessible with current computing power. 


Nonetheless, a number of codes to simulate global spherical convection are in use ( 

Miesch et al. 

2008; Ghizaru, Charbonneau &: Smolarkiewicz 20101 Kapyla et al. 2010, 2011; Hotta, Rempel &: 

Yokoyama 

2014 

). Most of these methods invoke the anelastic approximation (Gough 

1969 

), a 


low-Mach-number limit of the Navier-Stokes equation for stratihed media in which acoustic waves 
are hltered out. Removing acoustic waves substantially relaxes the Gourant restriction on the 


time step, thereby resulting in computational savings. Alternately, Hotta, Rempel & Yokoyama 


(2014) have developed a technique in which they artihcially reduce the sound speed, allowing weak 


compressibility to be incorporated (in contrast to the anelastic limit). 


3. Helioseismology of solar convection 


Helioseismology provides an independent and powerful means to constrain properties of interior 
convection. The Sun and its main-sequence cousins support a spectrum of oscillations, continuously 
excited by vigorous turbulent motion. The descending plumes in the intergranular lanes accelerate 
as they fall inward and are thought to generate shocks, thereby channeling their energy into acoustic 


radiation {e.g. Leibacher & Stein 1971). While the generation of acoustic waves is a nonlinear 


process, the subsequent wave propagation is well described by linear theory owing to the fact that 
solar acoustic mode amplitudes are substantially smaller than sound speed. Solar oscillations have 


resulted in accurate inferences of the structure of the Sun (Ghristensen-Dalsgaard et al. 1996), its 


internal rotation {e.g. Schou et al. 1998), and other large-scale flow features such as meridional 
circulation (Giles et al. 1997) and toroidal oscillations (Howe et al.||2000). 


Seismic inference involves measuring oscillation frequencies or wave travel times, and building 
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models of a medium in which the predicted properties of wave propagation would match observa¬ 
tions. Solar oscillations were discovered by Leighton, Noyes & Simon (1962), who speculated that 
these waves could be used to infer the properties of the surface layers. Subsequent measurements 
showed that the Sun exhibits acoustic waves that are trapped in the radial direction (for a review, 
see e.g. Christensen-Dalsgaard 2002), which could be used to probe the structure and dynam¬ 
ics of its deep interior. Contemporary helioseismology deals with the study of non-axisymmetric, 
three-dimensional perturbations due to magnetic helds and convective flow systems {e.g. Gizon, 


Birch Sz Spruit 2010). In subsequent discussions, we invoke the term “horizontal” which refers to 


the latitude-longitude directions and “radial” to the spherical radial direction. In many instances 
we perform the seismology of small regions, small in comparison to the solar radius. The formal 
interpretation of these measurements is expressed in plane-parallel Cartesian geometry, and the 
term “radial” is replaced by “vertical”. Two robust quantities that are typically inferred from seis¬ 
mology are the horizontal flow divergence and radial component of flow vorticity, i.e., V/^ • u and 
V/iXu, where V/i is the horizontal component of the gradient operator and v is the flow velocity 
(see Langfellner, Gizon & Birch 2014, for a description of how these measurements are made). 


3.1. Seismic measurements 


The Helioseismic and Magnetic Imager (HMI; Schou et al. 2012), an instrument onboard 
the Solar Dynamics Observatory launched by NASA in 2011, provides an important source of 
helioseismic observations. Among other measurements, HMI takes 4096 x 4096-pixel images of the 
line-of-sight velocity at the solar surface every 45 seconds (Doppler shifts of the Fe-I absorption 
line). Ground-based instruments such as the Global Oscillation Network Group (GONG) also 
provide streams of observations. Seismic waves are continually excited by turbulent convection in 
the Sun and consequently the wavefield is stochastic. Following Duvall et al. (1993), we measure 
the two-point correlation function C of the wavefield 4>{r, t) from finitely long observational records. 


Cafi{t) = ^j^ (j){ra,t')(j){rf},t'+ t)dt', (1) 

where T is the duration of the observation, t is the correlation time lag, and Tq, and are two 
locations on the surface. The Sun is continuously evolving with time and so averages over some finite 
temporal extent T imply that the dynamics of interest is quasi-static over T. For each pair of points, 
the travel time for waves propagating from to r^, denoted by Tap, can be extracted from the 


t > 0 part of the cross-correlation function ( 

Gizon Sz Birch 

2002 

2005 

). In turn, travel times may 

be formally interpreted (IGizon Sz Birch||2002 

Hanasoge et al.||20i: 

) in terms of products of Green’s 


functions of operator ([^ and an inverse problem can be posed. The travel time differences STap = 
Tap — Tpa are particularly sensitive to flows (as opposed to isotropic wave-speed perturbations) and 
are thus used as input to the inverse problem. 
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3.2. Seismic modeling 


Wave propagation in the Sun can be very effectively treated as a linear problem owing to 
their small amplitudes in comparison to the sound speed (wave amplitudes are on the order of 
cm/s compared to km/s sound speeds; Christensen-Dalsgaard 2002). The displacement vector ^ 
associated with small-amplitude waves satisfies {e.g. Lynden-Bell &: Ostriker||l967) 


p{dt - V • V)2^ - ^ ^ . Vp) - g V • (p^) = S, 


( 2 ) 


where p, p and c are the background density, pressure, and sound speed at position r and time t and 
g = —g{r) f is the acceleration due to gravity. In equation (§ , the first two terms model sound¬ 
wave propagation in a stratified medium and the third term is the effect of buoyancy. Advection 
by a background flow v is included. In order to facilitate seismic analyses, we assume that the 
background properties such as density, flows, pressure and sound speed vary on time scales much 
longer than the acoustic wave period and the observation duration. The source of excitation S 
represents random forcing by small-scale turbulence (granulation). 


Seismology is a technique to infer the wave speed distribution in the medium. In a conducting 
plasma, there are two basic types of wave speeds: an isotropic sound speed and magnetic Alfven 
velocity (e.g., Hanasoge et al. 2012). Additionally, flows in the medium advect acoustic waves in 
a manner that breaks the directional symmetry of propagation. Waves travel faster along the flow 
than against the flow. For instance, if the travel time between two points were to be r, the travel¬ 
time perturbation 6t introduced by a subsonic flow v is approximately given by 8t ~ —tv ■ k/c, 
where k is the unit wavevector. Thus a wave propagating along the direction of the flow will speed 
up, thereby reducing the travel time, and vice versa. 


Formally, seismic inference requires solving an inverse problem to estimate the coefficients of 
the partial differential equation Q, such that the difference between observation and prediction 
is minimized in some prescribed measure (such as the L 2 -norm difference). Given a model of the 
medium m = (p,c,v ■ ■ ■), we establish a cost functional y in terms of the travel times r, and thus 
which in turn is affected by the medium according to equation Q. An example of a travel-time 
misfit functional is 

( 3 ) 

a,P 

where the superscript ’obs’ denotes observations and Tap the predicted travel time between points 
Ya and r^. In this case, the goal would be to minimise the difference, measured in an L 2 sense, 
between the observed and predicted travel times. We consider the variations of x = x(^)j 

8x{m) = x(^ + ^'1^) — -|- 0(||5m|p). (4) 

To ensure convergence to the solution, i.e. 6x < 0, we choose 6m based on the negative gradient of 
X with respect to m, and use techniques such as steepest descent or nonlinear conjugate gradient 
to obtain improved estimates of m. Methods of optimization constrained by partial differential 
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equations (also widely used in aerodynamic shape optimization, e.g. 

Jameson 

1988 

) can be used 

to arrive at a best-fit seismic model of the medium ( 

Hanasoge et al. 

2011 

). In this process it is 


particularly important to take account of the correct noise matrix, which can be directly estimated 
from observations or using a model (Gizon &: Birch 2004 Fournier et al. 2014). More traditional 
approaches involve solving a linear inverse problem (the first iteration in the above problem); such 


approaches rely on careful choices for regularization terms {e.g. Svanda et al. 


2011 


Jackiewicz et al. 


2012 ). 


3.3. Weak convection on large scales 


An important first step in understanding convective flows is to provide estimates of their 


kinetic energy. To address this problem, Hanasoge, Duvall & Sreenivasan (2012) analyzed HMI 


observations of the seismic wavefield using techniques of time-distance helioseismology (Duvall 


et al. 1993[ ). These raw data were downsampled, cross correlations were measured according to 
equation 0- spatially averaged, and travel times measured. The interested reader may refer to 
Duvall| (|2003) and Hanasoge, Duvall &: DeRosa (2010) for details; the study by Hanasoge, 


e.g 


Duvall &: Sreenivasan (2012) inferred velocities at depth for spherical harmonic degrees ^ < 60 to 


be smaller than those predicted by MLT and numerical simulations (Miesch et al. 2009; Ghizaru, 


Charbonneau &: Smolarkiewicz||2010 Kapyla et al.|[20T0 ) by one to two orders of magnitude. They 
are also below the surface (Doppler) velocities {e.g. Hathaway et al. 2000; Rieutord et al.pOOS ), 
which are themselves below the simulated values. These results are shown in Figure Also 
shown in the figure are measurements of giant cell velocities on the order of 8-20 m/s inferred 


from supergranulation (SG) tracking by Hathaway, Upton &: Colegrove (2013). The corresponding 
kinetic energy is about an order in magnitude larger than seismic constraints at depth t/Rq = 0.96. 

Taken at face value, the seismic and giant-cell constraints of Figure [^suggest that our under¬ 
standing of thermal transport in the convection zone is incomplete. Because of the limited resolution 
of simulations and the weaker density gradients in them, the treatment of the highly asymmetric up 
and down flows is possibly not in the correct regime. Thermal transport of a solar luminosity can 
still be explained by decoupling the up and down flows and assuming that descending intergranular 


plumes make it all the way to the base of the convection zone (Spruit 1997) 


A related challenge lies in explaining the convective power fall-off at low spherical-harmonic 


degrees. Using phenomenological models. Lord et al. (2014) attempted to derive a set of conditions 


that would result in this fall off. In line with growing evidence supporting weak convective motions 
on large scales, they conclude: “we found that reducing the convective transport role of large-scale 
modes (by employing an artificial energy flux at all depths below 10 Mm which reduces the deep 
rms velocities by a factor of ~ 2.5) can significantly improve the match between the coherent struc¬ 
ture tracking spectra wof the simulations and observations. These separate lines of evidence all 
suggest that the Sun transports energy through the convection zone while maintaining very low 
amplitude large-scale motions. Something is missing from our current theoretical understanding 
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Fig. 5.— Upper bounds on kinetic energy E^j, of longitudinal velocities versus spherical-harmonic 
degree, i. We define at radius r such that (u^)/2 = X]£>o where the expectation value 

is approximated by a horizontal average. The grey region shows the helioseismology bound, based 
on T = 96-hour samples of HMI observations (Hanasoge, Duvall & Sreenivasan 2012). These 
constraints are several orders in magnitude smaller than numerical simulations of global convection 
(Anelastic Spherical Harmonic simulation (ASH); Miesch et al. 2009), suggesting that our current 
modeling of large-scale convection in the Sun is incomplete. The various curves denote convective 
energy spectra at different depths: seismology corresponds to v/Rq = 0.96, ASH to v/Rq = 
0.97, and Stagger simulations to v/Rq = 0.98. The horizontal black line is a theoretical lower 
bound based on global dynamics arguments Miesch et al. (2012), assuming mode equipartition 
over i < 750. The red curves are surface spectra based on HMI observations of granulation and 


supergranulation (SG) tracking. The SG tracking spectrum is based on data from Hathaway, Upton 


& Golegrove 

(2013), courtesy of David Hathaway. Adapted from 

Gizon &: Birch 

(2012 
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of solar convection below ~ 10 Mm”. Of particular relevance to models of convection are the 
matching of the amplitude of convective motions at all scales at the surface and accurately match¬ 
ing the ratios of power spectra at different depths. Both of these factors are strong tests of the 
mechanism of convection; the constraints in Figure suggest that the amplitude and depth-scaling 
of the convective spectrum are in question. More work is needed to have greater confidence in the 
helioseismic inferences on deep convective velocities. In this regard, several ideas for improvement 


have been identified by Gizon & Birch (2012). In particular, the treatment of the interaction of 


solar acoustic waves with turbulence is not fully accounted for in current analyses {e.g. Hanasoge, 


Gizon & Bal 2013). 


The ratio of the horizontal dimension to the radial dimension of the convection region is the 
aspect ratio. The effective aspect ratio of the convective zone in the Sun is large (of order 10) 
but not in most high-Rayleigh-number laboratory studies (of order 1). Nevertheless, some large- 
scale structural features of convection are shared between laboratory experiments and observations. 
Higher aspect ratios support large-scale motion that extends over the entire extent of the apparatus. 


but it is relatively weak; see Niemela & Sreenivasan (2006a). 


3.4. Flow divergence and vorticity 


Local helioseismology gives access to several useful quantities that describe turbulent con¬ 
vection in the near-surface layers: horizontal flow velocities, horizontal divergence, and two-point 


velocity correlations. Langfellner, Gizon & Birch (2014) measured the radial vorticity by com¬ 


puting the wave travel time around a closed contour. These quantities can help us characterize 
the dynamics and topology of rotating convection in the upper convection zone (from the surface 
down to a few Mm). For example, measurements of the radial vorticity provide a means to test 
models of the influence of rotation on convection {e.g. Egorov, Rudiger &: Zieglerj2004). Away from 


the equator, Langfellner, Gizon & Birch (2014) observe a clear correlation between radial vorticity 


and horizontal divergence. The latitudinal dependence of this correlation, proportional to HsinA, 
where H is the solar rotation rate at latitude A, is consistent with the action of the Goriolis force 
on convective flows. 

In addition to local helioseismology, surface flow velocities can be measured with the local cor¬ 
relation tracking (LGT) of granules. The idea is to use granules as tracers of larger-scale flows. The 
LGT technique is very useful as it provides an independent measurement of horizontal convective 
velocities at the solar surface. The method involves cross-correlating small patches containing a 
few granules in consecutive images measured in intensity (see Figure [^. The time lag between 
consecutive images (typically 45 s) is much less than the lifetime of granules (10 min), thereby 
avoiding contamination from granulation evolution. Gomparison of line-of-sight projected LGT 


velocity maps with Doppler and local helioseismology maps shows excellent agreement (Duvall & 


Gizon 2000 De Rosa, Duvall &: Toomre 2000). The spatial resolution of local helioseismology 


cannot be better than half the wavelength of seismic waves (> 5 Mm), while the resolution of the 



























- 17 - 


LCT technique of the order of few granules (> 2 Mm). Both methods are able to probe spherical 
harmonic degrees less than about 300. 


3.5. Influence of rotation on convection 


The effects of rotation on convection depend on the Rayleigh and Coriolis numbers. The 
highest Rayleigh numbers at which rotating convection has been studied in the laboratory are 
10^^ < Ra < 4.3 X 10^® (Niemela, Babuin Sz Sreenivasan 2010). The main result of this study 
is that modest rotation at high Rayleigh numbers has little effect on heat transport; if anything, 
rotation lowers the Nusselt number modestly. The vortical structure of the flows can however be 
affected quite significantly. 


Convection in the Sun is influenced by rotation in a complex way. Large effects of rotation 
might have resulted in standard Taylor-Proudman convective columns parallel to the rotation axis. 
In contrast, a rather small effect is observed. At smaller scales, the effects of rotation on convection 
become subtle as the Coriolis number becomes less than unity. The supergranulation flows are 
associated with a net kinetic helicity, Lfkin = (V A v) ■ v, and the horizontal divergence is anti¬ 


correlated with the radial vorticity in the northern hemisphere (Gizon & Duvall Jr. 2003; Gizon, 


Birch & Spruit] 

2010 

Langfellner, Gizon &: Birch 

2014 


and switches sign in the southern hemisphere. 


The positions of thousands of supergranules (in a latitude strip) can be read from the divergence 
maps using an image-segmentation algorithm. Averaging over outflow (or inflow) centers helps in 
lowering the noise level and in spatially resolving vorticity. Observations are shown in Figure 8 
for the LCT technique. In the north, outflows regions are associated with clockwise motion, while 
inflows are associated with anti-clockwise motion. The azimuthal velocity in the vortices is on the 
order of 15 m/s, i.e. 20 times smaller than the horizontal outflows. The patch of radial vorticity at 
the center of supergranules is more extended than in converging regions, where we observe a sharp 
localized peak in vorticity (width ~ 5 Mm). 


4. Future challenges 

4.1. Differential rotation and Reynolds stresses 

The outer convective envelope rotates differentially as depicted in Figure The inference 
of differential rotation in the convective envelope by helioseismology was remarkable because it 
revealed a profile that departs significantly from Taylor-Proudman balance (co-axial cylindrical 
iso-rotation). In the limit of weak Reynolds stresses, it is possible to invoke thermal-wind balance 
for a rotating spherical fluid, which balances the latitudinal entropy gradient and the resulting 
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rotation profile as 


9 d{S) 
dz Cp r cos 6 dO ’ 

where 0 is co-latitude, r is radius, z = rcosO is the distance along the rotation axis, = Cl{r,6) 
represents the rotational shear in the convective envelope, g is gravity, (5) = {S){r,6) is the 
temporally and longitudinally averaged entropy and Cp is the specific heat at constant pressure. 
For fluid spheres for which the entropy is a constant, we recover the standard Taylor-Proudman 
balance, in which the rotation rate H is purely a function of the distance from the axis and contours 
of constant rotation would be parallel to the rotation axis. However, it is seen from Figure|^that the 
observed rotation rate is constant along conical contours, i.e. dQ/dz / 0, suggesting the prevalence 
of a latitudinal entropy gradient. 


Equation ([^ by itself is unable to fully constrain the rotation rate since it allows for arbitrary 
terms that have no dependence on z to be added to The determining factor is the Reynolds- 
stress forcing, known in atmospheric-physics parlance by gyroscopic pumping (Miesch et al. 2012). 
Correlations between convective velocities, 


Rij = {Viv'j), 


( 6 ) 


yield the Reynolds-stress tensor TZij = pRij, where u' are the velocity fluctuations ((u') = 0) and 
the angular brackets indicate longitudinal averaging. Measuring the Reynolds-stress distribution 
in the Sun is thus a major challenge and one with important consequences. The weak motion of 
large-scale convection (see section 3.3) already suggests that we may expect some unknown issues 
to arise. 


Models of differential rotation were developed e.g., by Kitchatinov &: Rudiger (2005) and 


Rudiger, Kiiker &: Tereshin (2014), whereby angular momentum transport is due to small-scale 


rotating turbulence. In such mean-field models, transports angular momentum toward the 
equator as required to produce solar-like differential rotation. The model predicts that Rg^ is the 
sum of a viscous term (Kitchatinov, Pipin & Riidiger [1994 ) and the so-called A effect (anisotropy 
of turbulent convection caused by rotation, see Kichatinov &: Rudiger 1993), with a latitudinal 
dependence of the form Rg^ ~ C cos0sin^ 9, where 9 is co-latitude and C depends on the relevant 
Coriolis number. Hathaway, Upton &: Colegrove (2013) measured giant-cell velocity correlations 
at the surface with Cqc ~ 50 m^s“^, while seismic observations of supergranulation indicate a 
proportionality constant Csg 


— 10^ m^s ^ 


(see Figure 10 r). Although the relative effect of 
rotation on supergranulation is less than for giant cells (Coriolis number smaller), supergranulation 
velocities have larger amplitudes and thus measurable anisotropic Reynolds stresses. In the case of 
the (deep) giant cells, the A effect (positive C) dominates over the viscous term (negative C) and 
transports angular momentum towards the equator. Closer to the surface, however, the viscous term 
dominates in the supergranulation layer (Figure [T^). Observations and model are in remarkable 
agreement. 
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The model of Kichatinov & Rudiger (1993) not only predicts the components of the Reynolds- 
stress tensor, but more generally the two-point velocity correlations 


= Wiir'Wjir' + r)), 


(7) 


where r is the spatial lag and angular brackets denote an average over r' in a patch of the Sun around 
latitude A. The perturbation to Re(i,{r) introduced by rotation, denoted by 6Rg(i,{r), is expected to 
be mirror symmetric and to depend on the latitude through the local Coriolis number. In order 
to check the predictions of the model, measurements of the two-point correlation were obtained 
using time-distance helioseismology (Figure [TT|). We find that 6Rxy has the correct symmetries, 
has different signs in the northern and southern hemispheres, and vanishes at the equator. 


One of the interesting dynamical properties of the upper layers of the Sun is the existence of a 
strong radial gradient of rotation {e.g. Schou 2003). In this layer, which occupies about 5% of the 
solar convection zone, the rotational velocity increases with growing depth at the rate ~ 2 m/s per 
Mm. The latest measurements by Barekat, Schou & Gizon (2014) using HMI observations show 
that the logarithmic gradient of the rotation rate (dlnll/dlnr) near the surface is close to —1 from 
the equator up to 60° latitude. 


Both local heliseismology and LCT present serious possibilities for measuring Rij near the 
surface. Future work aims at including a larger range of spatial scales and probing deeper in the 
convection zone. The r-components of the Reynolds-stress tensor are more difficult to measure due 
to the fact that radial velocity fluctuations n/ are very small. 


4.2. Large-scale magnetism 


The Sun supports a global magnetic field that undergoes cyclic polarity reversal on a time-scale 
of roughly 11 years (in comparison, the geodynamo is erratic and reverses polarity on timescales 
of thousands of years). Elucidating the mechanism of the solar dynamo is a central goal of solar 


physics. A widely held description of the solar dynamo is the Babcock-Leighton mechanism (Dik 


pati Sz Charbonneau|||19^ ), which relies on differential rotation to intensify magnetic fields buried 
at the base of the convection zone. Rotational shear would strengthen these fields till they become 
magnetically buoyant and rise to the surface, where they would ostensibly form sunspots and other 
large-scale magnetic complexes. Magnetic tubes are much like rubber bands and the amplification 
process would require angular momentum to be extracted from rotational shear. However, dif¬ 


ferential rotation is observed to be remarkably constant (with small fluctuations, e.g. Howe et al. 


2000), and this would indicate that shear in angular momentum is being continuously replenished 


by convective Reynolds stresses, as is generally believed. 
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4.3. Super granulation 


Supergranulation is a prominently visible scale of convection on solar surface (Figure[^ though 
the thermal contrast is weak; they are only slightly warmer, likely less than 3 K at their centers {e.g. 


Meunier, Tkaczuk & Roudier 2007; Goldbaum et al. 2009). Consequently, the driving mechanism 


of supergranulation and its role in global solar dynamics is a matter of some debate (see Rieutord 


et al. 2010, for a review). 


It is not known definitively whether supergranules are deep or shallow structures (connected 
to deep or shallow convection processes). Local helioseismology shows that they penetrate at least 


down to a depth of 5 Mm {e.g. Sekii et al. 2007 Jackiewicz, Gizon & Birch 2008). However, one 


cannot exclude the likelihood that supergranulation downdrafts remain coherent over much greater 
depths. The seismic inference of the internal structure of supergranulation provides an effective 
means of discriminating between theories. Thus, supergranulation has been a subject of central 


interest in helioseismology {e.g. Kosovichev &; Duvall 1997| Duvall & Gizon 

[2000 

Braun, Birch & 

Lindsey 

2004 

Svanda et al.||2011 

Dombroski et al. 

2013 

)• 


Duvall, Birch &: Gizon (2006) pioneered the statistical study of ensembles of supergranules. 


based on the assumption of horizontal homogeneity. Analyzing and averaging wave travel times 


over thousands of supergranules, Birch et al. (2006) obtained precise estimates for the travel times 


associated with the average supergranule. These travel times can ostensibly be used to create 
accurate models of supergranules. However, seismic inferences have been unable to arrive at a 
unified model that successfully explains observed travel times. A number of these efforts are at 


odds with each other, and analyses by Duvall & Hanasoge (2013), Svanda (2012) and Duvall, 


Hanasoge &: Ghakraborty (2014) have suggested a previously unanticipated result, namely that 
supergranular flows peak at a depth of 2-3 Mm, reaching up to 800 ms“^ in horizontal flow speed. 


In contrast, Birch et al. (2006) and Woodard (2007) infer supergranular flows that decrease in 


magnitude with depth. Much work, in terms of refining solutions to both the forward and inverse 
problems in local helioseismology, is needed to satisfactorily address these issues. 


Duvall (1980) showed that the supergranulation pattern rotates 3-4% faster than the local 
plasma. This phenomenon is known as the superrotation of supergranules. The evolution of the 


supergranulation pattern was analyzed by [Gizon, Duvall Sz Schou (2003) using long sequences (two 
to three months each year) of seismic maps of horizontal divergence of the flows (see space-time slices 
in Figure [T^. Fourier analysis of these data revealed that the supergranulation pattern has wave¬ 
like properties with oscillation periods in the range of 6-9 days and excess power in the direction 


of rotation, in accordance with an independent analysis based on direct Doppler imaging (Schou 


2003). Gizon, Duvall &: Schou (2003) suggested that supergranulation is traveling-wave convection: 


the supergranulation pattern can be modeled by a superposition of modes (with anisotropic power) 
that travel at a phase speed of about 65 m/s. 


Since Chandrasekhar (1961), who investigated the linear stability of thermal convection, it is 


known that time-dependent oscillatory modes can develop in thermal convection in the presence 
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of either rotation or the magnetic held. Furthermore, the presence of shear signihcantly inhuences 
convectively unstable modes, which become traveling modes and propagate horizontally. Adam 


(1977) studied a convectively unstable layer bounded below and above by stable layers and found 
that the phase speed of the traveling convection modes can cover the full range imposed by the 
shear how. Busse (2003, 2007) showed that rotating convection cells will in general exhibit a drift 
in the direction of rotation and generate a mean how. We note that the solar near-surface shear 
how (~ 70 m/s) is comparable to the phase speed of the supergranular modes (65 m/s). Whether 
this is a coincidence or a physical effect should be investigated, possibly via numerical simulation. 


5. Summary points 

• Convection in the Sun operates in an extraordinary parameter regime, Re 10^^, Ra 

1020,Pr ~ 10“®. 

• Outstanding challenges in understanding the Sun lie in estimating the Reynolds stresses and 
establishing the mechanism governing the emergence of large-scale, huid circulations (such 
as differential rotation and meridional circulation) and the sustenance of its global magnetic 
held. 

• Traditional means for understanding these challenges have relied on mixing length theory and 
important numerical simulations (which, however, do not resolve the how held even remotely 
well). 

• The Sun’s interior is optically inaccessible but supports acoustic waves, well approximated 
by linear theory. 

• Seismic waves are advected by the underlying convecting medium, thereby shifting their travel 
times, which have been used to infer motions of the medium. 

• Helioseismology has been used to constrain the amplitudes of convective velocities at large 
scales (i < 60) to be two orders in magnitude lower than current theory and computation, 
throwing into question our understanding of thermal and angular momentum transport in 
the solar convection zone. 

• Seismology is also used to image intermediate scales of convection such as supergranulation 
(£ « 120) and to measure the anisotropic Reynolds stresses that control the global dynamics 
of the solar convection zone. 


6. Future issues 

• A complete seismic analysis of convective velocities that addresses a range of depths, from the 
mid-convection zone to the near-surface layers while simultaneously accounting for sources 
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of theoretical and observational errors is critical. Although a difficult problem, it will be a 
significant advance in placing seismic constraints on convection on a rigorous footing. 

• Numerical simulations of solar convection miss some critical ingredients since they are unable 
to accurately reproduce solar differential rotation; it is unlikely that one can push the simula¬ 
tions to the appropriate parameters, so one needs to have some independent and quantitative 
means for contrasting the numerical data. The path forward likely lies in the identification of 
important and relevant physical processes that control the dynamics of the convection zone 
and simulating these phenomena. 

• Laboratory experiments can lend significant insight towards appreciating the properties of 
highly stratified and non-Boussinesq convection, but studies of this sort are very sparse. 
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Annotated references 


Dnvall et al. (1993) introduced the method of time-distance helioseismology, widely used 


to image three-dimensional structures and flows in the solar interior. 


Hanasoge, Duvall & Sreenivasan (2012) bounded solar large-scale convective velocities. 


constraining them to be two orders in magnitude smaller than theory. 


Schou et al. (2012) described the purpose and workings of the Helioseismic and Magnetic 


Imager (HMI) onboard the Solar Dynamics Observatory (SDO), launched by NASA in 2010. 


Miesch (2005) gave an overview of numerical simulations of solar convection. 
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Gizon, Birch & Spruit (2010) summarized methods and results of local helioseismology. 


Niemela et al. (2000 ) performed convection experiments using liquid Helium at extremely 


high Rayleigh numbers. 


8. Related Resources 

• HMI data are publicly available at http://hmi.stanford.edu and from the German Data Center 
at http://www2.mps.mpg.de/projects/seismo/GDC-SDO 


9. Sidebar 


Alongside Section 1.1 

Much like the Sun, the Earth also possesses a convection zone (albeit in a vastly different parameter 
regime), termed the mantle, where overturning convection is driven by heating from the core. 
Surface cooling results in the formation of a crust and at subduction zones (such as the Pacific 
“ring of hre”), pieces of the crust break off and fall back into the interior. Terrestrial seismic 
imaging suggests that these pieces of crust (likened to the descending plumes in the Sun) reach the 
core-mantle boundary intact (e.g. Woodward &: MastersIlQ"^ ). Based on this picture, it is tempting 
to speculate that descending cool plumes formed at the solar surface maintain their identity on their 
journey to the base of the convection zone. This suggests that convection is not well mixed in the 
solar interior, indicating small convective velocities, in line with seismic and related inferences. 
It must, however, be recognized that the parameter regimes are starkly different; the Reynolds 
number is very small within Earth’s mantle and the Prandtl number is on the order of 10^ (and 
dramatically varying with the depth). 


10. Acronyms and definitions 

• ASH: Anelastic Spherical Harmonic code, that simulates global convection in the Sun in the 
anelastic limit. 

• LCT: Local Correlation Tracking is a technique that uses granules to trace larger-scale flows. 
It provides an independent measurement of horizontal convective velocities at the solar sur¬ 
face. 

• GONG: Global Oscillation Network Group, a network of ground-based telescopes around the 
world that observes the Sun continuously. 

• HMI: Helioseismic and Magnetic Imager, the helioseismology instrument onboard SDO 
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• SDO: Solar Dynamics Observatory, a NASA satellite in geosynchronous orbit launched in 

2010. 

• Rayleigh Number: Ra is the ratio of buoyancy to viscous forces multiplied by the ratio of 
viscous to thermal diffusivities. Beyond a critical value of the Rayleigh number, convection 
sets in; below this value, conduction is the primary mechanism of thermal transport. The 
mathematical definition can be more complicated for stratified fluids. 

• Reynolds Number: is the ratio of inertial to viscous force. Re = UL/v, where U is the char¬ 
acteristic convective velocity, L the characteristic length scale and v the kinematic viscosity. 

• The Coriolis number is the ratio of convective turnover time {L/U) to rotation period 27r/D, 
where D is the rotation rate, i.e. Co = LQ/(2 ttU). R is inversely related to the Rossby 
Number, which is the ratio of inertial to Coriolis force, Ro = U/{2LQ), by the relation 
Co ~ l/(47riio). 

• Prandtl Number: the ratio of viscous to thermal diffusivities, Pr = 
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Time-distance helioseismology Granulation tracking 



Fig. 6.— Measurement techniques in local helioseismology and local correlation tracking (LCT). 
In time-distance helioseismology, travel times of seismic waves measured between a central point 
and a concentric annulus inform us about the horizontal divergence of the flow, while travel times 
measured along a closed contour are related to the vertical vorticity of the flow. The motion of 
granules observed in time series of intensity images inform us about horizontal flows at the solar 
surface. 
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Fig. 7.— Horizontal divergence (top row) and radial vorticity (bottom row) of solar flows using 
HMI observations. The helioseismology (left panels, top 2 Mm) and LCT (right panels, surface) 
results are in good agreement: the correlation coefficients between the seismology and LCT maps 
are 0.93 for the divergence and 0.54 for the vorticity. The observation duration is eight hours and 
a low-pass filter was applied {I < 300). From Langfellner, Gizon & Birch (2014). 
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Fig. 8.— Vertical vorticity of the average supergranule in units of 10“® s“^ at latitudes —40°, 0°, 
and +40° (from left to right) from local correlation tracking analysis. The colors refer to the radial 
component of vorticity, the arrows to the horizontal flow components. The top row is for outflows. 


the bottom row for inflows. Adapted from Langfellner, Gizon &: Birch (2015). 
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Fig. 9.— (Left) Contour plot of solar interior rotation obtained from 19 years of data from 
the global network of ground-based observatories GONG. Black contours are at 10 nHz intervals. 
(Right) Rotational frequency as a function of fractional radius at selected latitudes. The thickness 
of the curves represents the formal (random) error. Note that due to resolution issues the tachocline 
(the shear layer at the bottom of the convection zone) is much thinner than shown here. Based on 
work by Howe (2009). Gourtesy of Rachel Howe. 
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Fig. 10.— {Left) Observed horizontal average Roff, = {v'gv'^) of co-latitudinal (southward) and 
longitudinal (prograde) fluctuating velocities as a function of solar latitude. The blue curve is for 
giant cells (local correlation tracking of supergranulation, Hathaway, Upton k, Colegrove||2013 ) and 
the data points and the red curve are for supergranulation (helioseismology, courtesy of Fournier, 
Gizon & Langfellner). The functional form cos0sin^0 is assumed. {Right) Rq^ at latitude 30° 


north as a function of fractional radius from a mean field model (cf. Riidiger, Kiiker & Tereshin 
2014). Units of m^s“^ are the same as in the left panel. Courtesy of Manfred Kiiker. 
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Fig. 11.— (Top panels) Two-point velocity correlation function, Rxy(r) = {Vx(r')Vy{r' + r)), as 
a function of horizontal lag r = {x,y), from local helioseismology and for latitudes —40°, 0°, and 
40° (from left to right). The coordinates x and y point in the directions of longitude and latitude 
respectively. The main features are due to non-rotating convection. (Bottom panels) Rotation- 
induced component of at latitudes —40°, 0°, and 40° (from left to right), obtained by taking 

the mirror-symmetric component of Rxy Based on the work by Fournier, Gizon, and Langfellner. 




























- 38 - 


12 


10 


^ 8 

cd 

^ 6 

a; 

6 
• 1=1 

^ 4 


2 


0 

0 10 20 30 40 50 

Carrington Longitude (deg) 



Fig. 12.— Evolution of the supergranulation pattern (divergence of horizontal flows) versus Car¬ 
rington longitude and time. The Carrington longitude is defined in a frame that co-rotates with 
the Sun. Each strip corresponds to a longitudinal section of the data centered at latitude 15° and 
1.68° wide. Supergranules (white regions of positive divergence) can be traced over several days. 
Several supergranules are long lived and are seen to travel faster (red arrow) than the local rotation 
rate. Erom Gizon (2006). 














